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We extend recent work by Tremblay, Turbiner, and Winternitz which analyzes an 
infinite family of solvable and integrable quantum systems in the plane, indexed by 
the positive parameter k. Key components of their analysis were to demonstrate 
that there are closed orbits in the corresponding classical system if k is rational, and 
for a number of examples there are generating quantum symmetries that are higher 
order differential operators than two. Indeed they conjectured that for a general 
class of potentials of this type, quantum constants of higher order should exist. We 
give credence to this conjecture by showing that for an even more general class of 
potentials in classical mechanics, there are higher order constants of the motion as 
polynomials in the momenta. Thus these systems are all super integrable. 
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I. INTRODUCTION 

Recently Tremblay, Turbiner and Winternitz [l| studied a family of quantum mechanical 
systems in two dimensions with Hamiltonian 

H = d r+-dr + ^ + a^ + ^\- + ^ Y -. (1) 

r r z r z cos z kv r A sin kv 

They showed that for k an integer this system is superintegrable, as is the corresponding 
classical analog. We will prove that if k is rational then the corresponding classical problem 
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with Hamiltonian 



H = p 2 r + \p 2 9 + ar 2 + b + ° (2) 



is also superintegrable. Indeed, it has two functionally independent constants of the motion 
that are polynomial in the momenta. One of these is always of second order corresponding 
to separation of variables in polar coordinates, viz 

L2= ^ + oo^ + sin^- (3) 

It is the other generating constant of the motion on which we will concentrate. 

To make this all explicit we consider the case k = 2 first and then give a general con- 
struction for rational k. Following this we consider a second potential 

{x + iy) k - 1 

where we also demonstrate superintegrability for k rational. These results strongly suggest 
that corresponding quantum systems are superintegrable. 

II. THE POTENTIAL OF TREMBLAY,TURBINER AND WINTERNITZ 

We look first at the k = 2 case. In Cartesian coordinates the classical Hamiltonian is 

H = A+ P l + a^ + ^)+b^±^ + c^l. (5) 
There are two fundamental constants of the motion: 

C\ = [xpy - yp x ) + 46— — + c- 



(x 2 — y 2 ) 2 x 2 y 2 ' 

C = (Pi - Vlf + [2ax 2 + 2b^±^- 2 - 2c^l]p 2 x 

+[ ~ Aaxy + m jx^y¥ p ^ + [2ay2 + 2b W^7T 2 + ^^^^ 

+a (x - y ) + — — + c — h 8ab— — + 2 



(x 2 — y 2 ) 2 x 4 y 4 r(x 2 — y 2 ) x 2 y 2 

The structure of the symmetry algebra is not worked out in [lj. We give the structure 
here. We set R = {Ci,^}. The Poisson algebra relations are 

{d, R} = 32(H 2 - 2C 2 )d - 64(6 + 2c)C 2 + 64(6 - c)H 2 - 128a6d - 128a6(6 + 2c) 
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{C 2 , R} = 32C 2 (C 2 - H 2 ) + \2SadH 2 - 384a 2 C 2 + 128a6C 2 - 64(6 + 4c)aH 2 
+256a 2 (2c - b)C x + 128a 2 (6 2 + 40c 2 + 206c). 
There is a Casimir constraint 

R 2 = QAC^iH 2 - C 2 ) - 646# 4 + 128(6 - c)C 2 H 2 - 64(6 + 2c)C 2 2 - 128aC\ 2 # 2 + 

256a 2 C 1 3 - 2hQabdC 2 + 128a(6 + 4c)if 2 C 1 + 256a 2 (6 - c)C 2 - 256a6(6 + 2c)C 2 

+256a(76c + 6 2 - 2c 2 )# 2 - 256a 2 (6 2 + 4c 2 + 206c)Ci - 256a 2 (2c + 6)(6 2 + 166c - 4c 2 ). 

From these relations we see that we have a closed Poisson algebra in the same sense as 
found for many well known superintegrable systems in two dimensions that are of second 

nnn 

order |2J, 13|, |4J] . In addition we can look for interesting one variable models of this algebra 
js|. One such model is obtained by choosing C% = c. If we do this then 

1 (4Cc - 4c 2 - C 2 + 166c) (4a6 + 9ac + AaC - E 2 ^ 2 



16 (C + 6 + 2c) 2 

Here (3 is the variable conjugate to c. 

For the quantum analogue of this system the Hamiltonian is 

H = 31 + d] + «(*" + y>) + bg±£L + 

with quantum symmetries 

d = (xd y - yd x ) 2 + 46 . f V \. 2 + c {x4 +f\ 



exp(- v / -C-6-2c/5). 



C2 = {d 2 x - d 2 ) 2 + (2ax 2 + 2b ^ + V 2 - 2c^-^l)d 2 x + 



[x* — y^y x*y z 



, , 86x7/ x _ _ , n 9 nl {x 2 + y 2 ) n (x 2 - u 2 ) . o9 
{-Aaxy + — — y -^)d x d y + (2m/ 2 + 26-^ ^ + 2c v - | ; )c> 2 



(x 2 — y 2 ) 2 (x 2 — y 2 ) 2 x 2 y 

+ (2ax - %)d x + (2ay - %)d y + a\x 2 - y 2 ) 2 + 2 * + ^ ~ f ] \ 
x 6 y A [x* — y ) x 4 y 4 

, x 2 y 2 26c ,1 L 

For these quantum operators there is corresponding closure given by the formulas ({•,•}, 
{•,-,•} are operator symmetrizers) 

R — [Ci, C 2 ], 



[d, R] = 32dH 2 - 32{C h d} + 64(6 - c + 2)H 2 - 64(6 + 2c + 4)C 2 - 128a(6 + l)C x 

-128a(6 2 + 26c + 46 + 6c + 4), 
[C 2 , R] = 32C\ - 32H 2 C 2 + 128adH 2 + 128a(6 + 1)C 2 - 64a(6 + 4c + 6)H 2 - 384a 2 C 2 
-256a 2 (6 - 2c - U)d + 128a 2 (-8 + 8c + 186 + 206c + 6 2 + 4c 2 ). 
There is also the Casimir operator 

R 2 = 32H 2 {d, d} - ^{d, d, d} ~ (64/3 + 128 7 + ^)C 2 2 + (128(0 — 7) + ^)H 2 d 



-(192 + 64/3)# 4 - 128aH 2 Cf - 128a(/3 + l){d, d} 
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128 



-a 



(12 7 + 3/3 + 50)# 2 Ci + 256aCf 



256 256 
•— a{U(3 + 44 + I87 + 3/3 2 + 6/37)C 2 - 256a 2 (2 7 - /346)C 1 2 + — a(42 + 22 7 + 40/3 
o o 

25fi 

+3/3 2 - 6 7 2 + 2\fii)H 2 - — a 2 (1527 - 88 + 182/3 + 3/3 2 + 12 7 2 + W(3^)d 



+ 



256 



(28O7 2 - 80/3 2 + 24 7 3 + 3207 + 48/3 - 4/37 - 84/37 3 - 54 7 /3 2 - 3/3 2 + 28) 



We now prove our central result that the classical Hamiltonian has, in addition to the 
obvious second order constant of the motion, another independent constant of the motion 
that is polynomial in the momenta. We use the results of a previous paper 6j. A similar 
approach was used by Verrier and Evans [7]. For the general potential we have in polar 
coordinates 



V = ar 2 + 



+ 
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r 2 cos 2 (k8) r 2 cos 2 (k$) 
In terms of the new variable r = e R the Hamiltonian assumes the form 



H = e- 2R (p 2 R + p 2 + ae 4K + 



4R 
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cos 2 (^) sin 2 (^) 



Applying the method of j^] to find the extra invariants we first need to construct a function 
M(R,p R ) which satisfies {M,H} = e~ 2R , or 



Aae m + 2He 2R )d PR M + 2p R d R M = 1. 



This equation has a solution 



M 



-.B 



where 



sinh£> = z- — . 2 „ ==)-, coshi? 
^H 2 - AaL 2 



2^L 2 e~ 2R p R 
^H 2 - 4aL 2 ' 
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and 

L 2 = pi + 7T7T~^\ + 



2 | g | 7 

2 ^ cos 2 (^) sin 2 (£;fl) 



and we also have the relation (which we can use to consider M as a function of R alone): 

p 2 R + L 2 + ae m - e 2R H = 0. 
We now need to find the corresponding function N{6,pq) which satisfies {iV, H} = e~ 2R , or 

+ _. ZmJ d Pe N ~ 2 Ped e N = 1 



v cos 2 (fc#) sin 2 (fc#)' 

where the prime denotes differentiation with respect to 9. This equation has a solution 
iV = - A where 

4VL2K 



'(L 2 -!3- if - 4/3 7 s](Li - 3 - if - 4/5 7 

The constant of the motion is M — N, and, since it is constructed such that {M — N, L 2 } 0, 
it is functionally independent of L 2 , jf|. 

From these expressions for M and N we see that if k is rational, k = E ( where p, q are 
relatively prime integers) then 



sinh(- Aip ^ L 2 [N — M)) — - smh(qA + pB), cosh(-4ip^ L 2 [N — M)) — cosh(qA + pB), 

each give rise to a classical constant of the motion which is polynomial in the momenta. 

This can be seen by observing that these constants of the motion can be expressed as 
factor functions of L 2 and H times a factor which is polynomial in the canonical momenta, 
via the relations 

(cosh x ± sinh x) n = cosh nx ± sinh nx, cosh(x + y) = cosh x cosh y + sinh x sinh y, 
sinh(x + y) = cosh x sinh y + sinh x cosh y. 

In particular, 



[n/2] 

coshnx = ^2 sinh 2 - 7 a; cosh n_2j x, 



n 



n 



[(n=l)/2] / 

sinhnx = sinhx 

i=i V 2j - 1 



sinh 2 - 7 2 x cosh™ 2j 1 x. 
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Thus if p, q are both odd then 

C 

cosh(gA + pB) 
sinh(gA + pB) 



{L 2 -(3- 7 ) 2 - A^}iym-AaL 2 ]p ' 
\fU D 



[.J(L 2 - - 7)2 - 4 / 3 7 ]<V# 2 - 4aL 2 }p ' 
where C, D, are polynomial constants of the motion of orders 2(p+q), 2(p+q) — 1, respectively. 
If one of p, q is odd and the other even, then 

VT 2 C 



cosh(gA + pB) 
sinh(gA + pB) 



[^{L 2 - (3 - 7 ) 2 - A^nVH* - 4aL 2 f 

jy 



[y/fa - (3 - 7 ) 2 - 4f3i]«WH*-4aL 2 ]p 
where C, D', are polynomial constants of the motion of orders 2{p + q) — 1, 2{p + q), respec- 
tively. ( We have in fact produced two extra constants of the motion whose degree differs 
by 1. This is easily understood by realizing that we have one extra constant and its Poisson 
bracket with L 2 .) For example, if p — 1, q = 2 we have (with C = —7 + (3 — L 2 cos 9), 

{e- 2R p R {L 2 sin 2 9p 2 e - C 2 ) - sm9(2L 2 e~ 2R - H)C 



2^[L 2 

cosh(2A + B) = 



sinh(2A + B) = 



[(La - (3 - 7 ) 2 - 4/? 7 ]V# 2 -4aL 2 
(2L 2 e- 2i? - #)(L 2 sin 2 0pg - £ 2 ) + 4L 2 sin 9e~ 2R p 6 p R C 



[(L 2 - /? - 7 ) 2 - 4^7] VL 72 - 4aL 2 
The bracketed quantities in the numerators are 5th and 6th order constants of the motion, 
respectively. 

III. A NEW POTENTIAL 

In addition to the potentials of Turbiner et. al. there is another family based on the same 
principle. To illustrate the properties of this family consider 

tt 2 |2| {x + iyf , , 

H = P x +P y + a (6) 

This Hamiltonian admits three constants of the motion: 

Ki = (p x - ipyf - j^—y[-(iy + 3x)p x + {-ix + 3y)p y ], 
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K 2 = {xp y - yp x )(p x - ipyf + ^— - °— p[(3a; 2 + 3ixy - 2y 2 )p 2 x - (2x 2 + 3ixy - 3y 2 )p 2 y 

i(x + 3iy)(iy + 3x)p x p y - a 



2 (x + iy) c 



(x — iy) 6 ' 

. , 2 y(3x 2 - y 2 ) 

K 3 = {xpy - yp x ) + 2ia- 



(x — iy) 3 
The Poisson algebra relations are 

{K u K 2 } = 3iK 2 , {K u K 3 } = &K 2 , 

{K 2 ,K 3 } = QtK 1 (K 3 + a), 

together with the constraint 

K\K 3 - K\ + a(K 2 - H 3 ) = 0. 

There is also a sixth order symmetry K\. There are a number of one variable models to 
consider for this Poisson algebra which help with the formulation of corresponding quantum 
problems viz. 



(1) : = c, ifi = - J- cos(6 v / cTa ( 5), K 2 = -iVaE 3 fi\n{<o^f^Ta(3) . 

V a + c 

aE 3 

(2) : K x = c, K 2 = 3ic 2 f3, K 3 = -8c 2 f3 2 + — - a. 

cr 

(3) : K 2 = c, K x = ^, K 3 = -9(c 2 + aE 3 )p 2 - a. 

We see that (1) indicates a realization of the quantum operators in terms of difference 
operators and (2) and (3) a realization in terms of differential operators. 
Proceeding to the quantum analogue we obtain the operators 

Ki = (d x - idyf + _ Q ^ 3 [~(^ + 3x)d x + (3iy + x)d y ], 

K 2 = (xd y - yd x )(d x - id y f + a . \i(2y 2 - 3ixy - 3x 2 )d 2 x - {3iy + x)(iy + 3x)d x d y + 

[x — iy) 6 

i(2x 2 + 3ixy - 3y 2 )d 2 - 2i(3iy + x)d x - 2{iy + 3x)d y - 8i] + ia 2 ^ "'' ' 



(x — iy) 6 ' 



K 3 = (xdy - yd x ) 2 + 2tay( !r ^ 



[x-iyf 



H = d 2 + d 2 + a- {x + w)e 



' X ' ~y 



V + 2/ 2 ) 4 ' 
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with the commutation relations 

[K x , K 2 ] = 3iK 2 , [K u K 3 ] = 6iK 2 - 9K X , 

[K 2 , K 3 ] = 3i{K 1 , K 2 } + *(27 + Qa)K x + 9K 2 , 

and the analogue of the constraint 

1 9 63 

~{K 1: K u K 3 } - 3K\ - i-{Kx, K 2 } + ( — + 3a)Kl - 3aH 3 = 0. 

A one dimensional model of this algebra is 

i d 
K 1 = , K 2 = —. 

3x ax 

d 2 d 
K 3 = -9x 2 — - 27x— - (9 + a + 9aE 3 x 2 ). 

We now look at the question of what the constants of the motion might be for more 
general potentials of type §6§. We consider the potentials 



(x + iy) k J 
(x - iy) f 



^ T oy 

v = a T-—rm- ( 7 ) 



As in the previous example it is convenient to pass to variables R and 6. In these 
coordinates the Hamiltonian and the obvious constant of the motion L assume the form 

tt _ (Pr + Pe + ae 2lke ) T 2 n jiiW T/ n 2ike-2R 
n — , L = p 9 + ae , V = ae 

Using the the usual prescription for obtaining the extra constant we need to look for 
solutions of 

2He 2R d PR M + 2p R d R M = 1. 

and 

-2tae 2ikd d P0 N + 2p e d e N = 1 

The new constant is then M — N. If k is an integer it is convenient to consider the solution 
such that -ik\fL(M - N) = A + kB where 

sinh A = ^=e~ lk \ cosh A = J^e~ ike , 
'a V a 



sinhjB = ^e" H , coshfi = J^e~ R . 
y/H V H 
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If k = 2 is rational then we consider sinh(qA + pB) and cosh(gA + pS) in order to obtain 
extra constants of the classical motion. For example, in the special case k = 2 we obtain 

sinh(A + 2B) = -^(2(po + i PR )L - p He 2R )e- 2R e~ 2W , 

coshM + 2B) = -^(2L - He 2R + 2%p e p R )e- 2R e- 2l \ 

where we need only consider these hyperbolic functions multiplied by y/aH in the first case 
and y/aH/ \J~L in the second, to obtain the polynomial solutions we seek. 



IV. CONCLUSION 

We have shown that for k rational all the classical mechanical systems ([I]) admit one 
second order constant of the motion as well as two others of higher order as polynomials in the 
momenta. This proves superintegrability and supports recent studies by TremblayTurbiner 



and Winternitz 



U, [8| of the potentials with k rational where it has been demonstrated that 
all the orbits are closed. We also studied a new class of systems (jZJ) and showed that again 
the systems are superintegrable and demonstrated how to find a maximal set of constants 
polynomial in the momenta. We provided some information about the structure of the 
symmetry algebras associated with all these systems. 
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